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Impacts 

• Coronavirus (CoV) was detected in 57 individual and 11 ground faecal 
samples from nine bat species in Taiwan, including severe acute respira¬ 
tory-related CoV, Scotophilus bat CoV 512 and Miniopterus bat CoV 1A. 

• Significantly higher detection rates of coronaviral RNA were found in 
female bats and Scotophilus kuhlii roosting in palm trees. 

• High nucleotide identities were shared between bat coronaviruses detected 
from the same bat species in Taiwan, China, and Philippines indicated the 
endemic circulation of bat CoV in local bat population by migratory bat 
species. 
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Summary 

Bats have been demonstrated to be natural reservoirs of severe acute respiratory 
syndrome coronavirus (SARS CoV) and Middle East respiratory syndrome 
(MERS) CoV. Faecal samples from 248 individuals of 20 bat species were tested 
for partial RNA-dependent RNA polymerase gene of CoV and 57 faecal samples 
from eight bat species were tested positive. The highest detection rate of 44% for 
Scotophilus kuhlii , followed by 30% for Rhinolophus monoceros. Significantly 
higher detection rates of coronaviral RNA were found in female bats and 
Scotophilus kuhlii roosting in palm trees. Phylogenetic analysis classified the posi¬ 
tive samples into SARS-related (SARSr) CoV, Scotophilus bat CoV 512 close to 
those from China and Philippines, and Miniopterus bat CoV lA-related lineages. 
Coronaviral RNA was also detected in bat guano from Scotophilus kuhlii and 
Myotis formosus flavus on the ground and had potential risk for human exposure. 
Diverse bat CoV with zoonotic potential could be introduced by migratory bats 
and maintained in the endemic bat population in Taiwan. 


Introduction 

In May 2015, Middle East respiratory syndrome coron¬ 
avirus (MERS CoV) caused 185 cases and 36 deaths in the 
Republic of Korea (Lee and Wong, 2015) when the virus is 
circulating mainly in Saudi Arabia since September in 2012 
(Zaki et al., 2012). One decade ago, severe acute respiratory 
syndrome (SARS)-CoV caused total 8096 cases worldwide 
and 346 cases with 37 deaths in Taiwan (Ksiazek et al., 
2003). Phylogenetic analysis showed that MERS-CoV clus¬ 
tered in the same subgroup with Tylonycteris bat CoV 
HKU4 and Pipistrellus bat CoV HKU5 (Lau et al., 2013). 


Later, bat CoV HKU4 was found capable of using the same 
cell receptor of MERS CoV (Yang et al., 2014b). MERS- 
related CoV was also detected in Nycteris bats of Ghana and 
Pipistrellus bats of Europe, Neoromicia zuluensis bats of 
South Africa, and Vespertilio superans bats of China (Annan 
et al., 2013; Ithete et al., 2013; Yang et al., 2014a). Rhinolo¬ 
phus bats have been suggested as natural reservoirs of 
SARS-CoV because the SARS related (SARSr) CoV found 
in Rhinolophus sinicus of Yunnan in China had 99.9% of 
sequence identity with human SARS-CoV and were able to 
use the angiotensin-converting enzyme 2 (ACE2) of 
humans, civets and bats as receptors to infect cells from 
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human, civets and bats (Ge et al., 2013). The outbreaks of 
MERS CoV re-emphasized the importance of bat-surveil¬ 
lance programs on CoV since the emergence of SARS CoV. 

The organization of CoV genome is in the order of 
5-polymerase-spike (S)-envelope (E)-membrane (M)- 
nucleocapsid (N)-polyA-3. The polymerase gene encodes 15 
to 16 non-structural proteins including RNA-dependent 
RNA polymerase (RdRp) gene, used commonly for phyloge¬ 
netic analysis (Woo et al., 2012). As the largest enveloped 
RNA viruses with single plus-stranded RNA genomes of 26 
to 32 kb, CoV are classified into four genera from Alpha- to 
Deltacoronavirus. The genus Betacoronavirus had four sub¬ 
groups from 2a to 2d, including SARSr CoV from human, 
civet cats, raccoon dogs and bats belonging to subgroup 2b 
and MERS CoV belonging to subgroup 2c (Agnihothram 
et al., 2014). Both Alpha- and Betacoronavirus have been 
found in the countries around Taiwan, including China, 
Philippines, Thailand and Japan. The most common bat spe¬ 
cies carrying bat CoV are Rhinolophus , Pipistrellus and Min¬ 
iop terus bats (Tang et al., 2006; Watanabe et al., 2010; 
Gouilh et al., 2011; Shirato et al., 2012), which distribute 
widely in Asia, including Taiwan. They could carry and 
transmit bat CoV into local bat population of Taiwan via 
migration even though little is known about the migratory 
patterns of bats. Closely related CoV can be detected in the 
same bat species living at locations separated by thousands of 
miles (Drexler et al., 2010) and different CoV species or gen¬ 
era can be found in different bat species living at the same 
roosting sites (Chu et al., 2006; Gloza-Rausch et al., 2008; 
Tong et al., 2009; Drexler et al., 2010). From a host switch, 
certain CoV species can be detected in totally different hosts 
(Gloza-Rausch et al., 2008). In addition, one bat species can 
harbour more than one CoV species (Chu et al., 2006; Tong 
et al., 2009). There are 11 endemic species and six endemic 
subspecies in 36 bat species in Taiwan (Wu et al., 2012; 
Csorba et al., 2014; Ruedi et al., 2015). Between two biodi¬ 
versity hotspots, Southeastern China and Philippines, Taiwan 
has an impressive biodiversity on unit area, 4.6 and 9 times 
higher than Thailand and Japan, respectively (Csorba et al., 
2014). It is critical for public health to understand the corre¬ 
lation between the prevalence of CoV and the dynamics of 
endemic and non-endemic bat species in Taiwan. 

Factors correlating with CoV infection include bat spe¬ 
cies, roosting location, habitat types, age, gender, physical 
condition and reproductive status (Gloza-Rausch et al., 
2008; Lau et al., 2010; Osborne et al., 2011). Because of the 
ability of fly and migrate and the tendency to live in the 
large sizes of social groups, bats are predisposed to main¬ 
tain and transport viruses to other mammals for the possi¬ 
ble outbreaks of epidemics (Bennett, 2006). Habitat types 
of bats include caves, trees and human-made constructions. 
The accumulation of bat faeces in schools or houses posts a 
potential risk for public health (Dobson, 2005). As a poten¬ 


tial zoonotic pathogen, CoV is regarded as particularly sig¬ 
nificant because its high mutation rate, high frequency of 
recombination, faecal-oral transmission, and environment 
resistance of faeces-immersed virus particles (Casanova 
et al., 2010; Graham and Baric, 2010; Geller et al., 2012). 
The investigation of RNA or infectious particles of CoV in 
the faeces accumulated on the ground under the roosting 
locations near or within human community, like Scotophilus 
bats in schools and houses, can be used to evaluate the expo¬ 
sure risk of CoV to human. Previous studies have shown that 
young and lactating female bats in maternity roosting sites 
had significant higher detection rate for bat CoV RNA 
(Gloza-Rausch et al., 2008; Lau et al., 2010; Osborne et al., 
2011). It is important to analyse the effects of bat physical 
properties on CoV infection for the risk assessment. 

The current study was to determine the prevalence of 
CoV in faecal samples collected from 21 bat species at 10 
sampling locations in Taiwan. Phylogenetical relationships 
of detected CoV were constructed by RdRp gene fragment. 
Factors contributing to CoV infection were analysed. The 
surveillance for CoV on bat reservoirs is important for con¬ 
trolling the potential threat of emerging CoV outbreaks in 
Taiwan and nearby areas. 

Materials and Methods 

Sample collection 

The protocol for capturing and sampling of bats was 
approved by the Chung Yuan Christian University Animal 
Care and Use committee and the agriculture bureau of regio¬ 
nal governments. All bat faecal samples were collected from 
August 2013 to October 2014 at 10 locations marked in 
Fig. 1 from the woods in Guihou (location 1, 25°11 / 24"N, 
121°40'41"E), Wulai (location 2, 24°5T50"N, 121°33'05"E), 
Dong’ao (location 3, 24°31 / 55"N, 121°49'53"E), Nan’ao (lo¬ 
cation 4, 24°3T9"N, 121°46 / 2"E), and Miaoli (location 5, 
24°25 / 35"N, 121°00 / 45"E), from the palm trees in Jhutang 
(location 6, 23°50 / 5T'N, 120°23 / 12"E) and Beigang (location 
9, 23°34 / 05"N, 120° 17'51 f 'E), from Ficus trees in Shuilin (lo¬ 
cation 8, 23°37 / 14"N, 120°15 / 40"E), and from the irrigation 
culverts in Dili (location 7, 23°48 / 24"N, 120°54'51"E) and 
Dongshan (location 10, 23°19 / 04"N, 120°25 / 27"E). A harp 
trap, butterfly net or mist net was used to catch bats. Species, 
body weight, gender, reproductive status and maturation sta¬ 
tus of bats were determined by experienced bat biologists at 
the Taiwan Endemic Species Research Institute. According to 
the swelling of nipples, female bats can be classified into lac¬ 
tating, weaning and non-lactating females. The mating status 
of male bats include mating, prepared to mate and non-mat¬ 
ing based on the production of semen by the descent testicles 
into scrotum. The bats with insufficient ephiphyseal-diaphy- 
seal closure in the metacarpal-phalangeal joints were consid¬ 
ered younger than 1 year old and others were considered 
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Fig. 1 . Map of East Asia (left) with enlarged view of Taiwan Island (right). The dots labeled with Arabic number indicated sampling locations 1 to 10. 
Name/habitat type/geographic coordinate: location 1/Guihou/woods/25°11 / 24"N/121°40 / 41 "E, location 2A/Vulai/woods/24°51 / 50"N/121°33'05"E, 
location 3/Dong'ao/woods/24°31 / 55"N/121°49 / 53"E, location 4/Nan'ao/woods/24°31 / 9"N/121 °46'2"E, location 5/Miaoli/woods/24°25'35"N/ 
121°00'45"E, location 6/Jhutang/palm trees in an elementary school/23°50'51 "N/120°23 / 12"E, location 7/Dili/irrigation culvert/23°48 / 24"N/ 
120°54'51"E, location 8/Shuilin/Ficus trees in an elementary school/23°37 / 14"N/120°15'40"E, location 9/Beigang/palm trees in a sugar factory/ 
23°34 / 05"N/120°17 / 51 "E, and location 10/Dongshan/irrigation culvert/23°19 / 04"N/120°25 / 27"E. 


adults. The captured bats were individually placed into cot¬ 
ton bags and allowed to produce faecal pellets. To clarify 
whether CoV RNA can be maintained in the surrounding 
environment, the accumulated faeces on clean plastic sheets 
overnight from Myotis formosus flavus at location 8 and those 
from Scotophilus kuhlii at location 6 and 9 were collected as 
ground faecal samples. Faecal samples were collected by 70% 
ethanol rinsed tweezers into tubes containing RNAlater® 
RNA stabilization solution (Qiagen, Flilden, Germany) for 
sample storage and transportation. 

Coronaviral sequence detection and analysis 

Viral RNA was extracted from the 140 -aL supernatant of 
each faecal sample by using the QIAamp Viral RNA Mini 
Kit (Qiagen) and reverse-transcribed into cDNA by using 
iScript Select DNA Synthesis kit (Bio-Rad, Hercules, CA, 
USA) with the mixture of random primers and oligo dT 
primers. Nested reverse transcription polymerase chain 
reaction (RT-PCR) with primers targeting a 440-bp RdRp 
fragment were performed on all samples according to previ¬ 
ous study (Poon et al., 2005). Sequences of the nested PCR 
products were determined by Sanger sequencing method 
(Genomics Company, New Taipei City, Taiwan). The con¬ 
sensus RdRp sequence (375-bp) of 24 bat CoV samples 


were aligned using the Clustal W method, and the phyloge¬ 
netic trees based on nucleotide sequences were constructed 
using the neighbour-joining model with a 1000 bootstrap 
test and pairwise deletion option in the maximum likeli¬ 
hood method through the mega6 program (www.megasoft- 
ware.net). The pairwise distance of nucleotide sequences 
was analysed using mega6 with the parameters of a 1000 
bootstrap test, a p-distance model, and pairwise deletion. 
Representative CoV sequences selected from NCBI for phy¬ 
logenetic analysis were listed in Table SI. 

Statistical analysis 

Comparisons of coronaviral RdRp gene detection results 
between groups were performed using Fisher’s exact test 
for categorical variables (bat species, gender, maturation 
status, reproductive status, sampling location and habitat 
type) and anova test for continuous variables (body weight) 
in the R project (R Core Team, 2013). P values of <0.05 
were regarded as statistically significance. 

Results 

Detection of bat CoV RNA in faecal samples 
A total of 248 individual and 64 ground faecal samples 
from 21 bat species were collected. Information regarding 
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Table 1 . Overview of bat faecal samples collected in Taiwan for detection of coronavirus 


Bat Species (abbreviation) 

Total No. (positive) 

Females (positive) 

Location 3 

Time (month/year) 

Family Hipposideridae 

Coelops frithii formosanus (CF) C 

1 (0) 

1 (0) 

5 

9/14 

Hipposideros armiger terasensis (H) b 

7(1) 

5(1) 

5, 10 

11/2013; 7, 10/2014 

Family Rhinolophidae 

Rhinolophus monoceros (R) b 

50(15) 

34(13) 

5, 7, 10 

11/2013; 7, 10/2014 

Family Miniopteridae 

Miniopterus fuliginosus (M) 

54(11) 

24(7) 

1, 5, 7, 10 

8, 9, 11/2013; 1, 7, 10/2014 

Family Vespertilionidae 

Barbastella darjelingesis (BL) 

2(0) 

0 

5 

7, 8/2014 

Harpiola isod on (HI) 

1 (0) 

1 (0) 

5 

8/14 

Kerivoula titania (K) 

10(1) 

6(1) 

3,4, 5 

9/2013; 7, 9, 10/2014 

Murina gracilis (MG) b 

2(0) 

0 

5 

9, 10/2014 

Murina puta (MP) b 

17(0) 

2(0) 

3,4, 5 

4, 7, 8, 9, 10/2014 

Murina recondita (MR) b 

4(1) 

3(1) 

4,5 

4, 6, 10/2014 

Myotis fimbriatus taiwanensis (MT) C 

22 (5) 

4(1) 

2, 10 

3, 10/2014 

Myotis laniger (MS) 

2(0) 

1 (0) 

5 

8/14 

Myotis rufoniger (MW) 

2(0) 

1 (0) 

3 

8/14 

Myotis secundus (MY) b 

9(0) 

6(0) 

4,5 

4, 6, 7/2014 

Pipistrel! us a bra m us (P) 

1 (0) 

1 (0) 

6 

9/13 

Pipistrellus montanus (PM) 

3(0) 

0 

5 

7/14 

Pipistrellus taiwanesis (PS) 

3(0) 

2(0) 

5 

8, 9/2014 

Plecotus taivanus (PT) b 

7(1) 

2(0) 

5 

7, 9/2014 

Scotophilus kuhlii (S) 

48 (21) 

28(11) 

6,9 

9/2013; 9/2014 

Submyotodon latirostris (ML) b 

Total individual faecal samples 

3(1) 

248 (57) 

0 

121 (35) 

4,5 

6, 9, 10/2014 

Scotophilus kuhlii (S) 

34(8) 

N.D. 

6, 9 

9/2013; 9/2014 

Myotis formosus flavus (MF) C 

Total ground faecal samples 

30(3) 

61 (11) 

N.D. 

8 

6/14 


a Name/habitat type/geographic coordinate: location 1/Guihou/woods/25°11 / 24"N/121°40'41 "E, location 2/Wulai/woods/24°51 / 50"N/121 °33'05"E, 
location 3/Dong'ao/woods/24 0 31 / 55"N/121°49 / 53"E, location 4/Nan'ao/woods/24 0 31 / 9"N/121°46 / 2"E, location 5/Miaoli/woods/24°25 / 35"N/ 
121°00'45"E, location 6/Jhutang/palm trees in an elementary school/23°50'51 "N/120 o 23'12"E, location 7/Dili/irrigation culvert/23°48'24"N/ 
120°54'51"E, location 8/Shuilin/F/a/s trees in an elementary school/23°37 / 14"N/120°15'40'T, location 9/Beigang/palm trees in a sugar factory/ 
23°34'05"N/120°17'51 "E, and location 10/Dongshan/irrigation culvert/23°19'04"N/120°25'27"E. 
b Taiwan endemic species. 
c Taiwan endemic subspecies. 

N.D. means 'not determine'. 


bat samples is listed in Table 1. Overall prevalence of CoV 
in all individual bats was 23% with detection rates of 10% 
in Kerivoula titania (K) to 43.8% in Scotophilus kuhlii (S). 
Coronaviral RNA can be detected in 17% (11/64) of 
ground faecal samples: 25% (2/8) from Scotophilus kuhlii at 
location 6 in September 2013, 23% (6/26) from Scotophilus 
kuhlii at location 9 in September 2014, and 10% (3/30) in 
Myotis formosus flavus at location 8 in June 2014. Of nine 
bat species tested positive for CoV RNA, five were endemic 
species and one was endemic subspecies in Taiwan 
(Table 1). The bats roosting in palm trees had significantly 
higher detection rates than those in irrigation culvert 
(P = 0.0243), followed by those in woods (P = 0.0127). Of 
10 bat species collected from more than two sampling loca¬ 
tions, Myotis fimbriatus taiwanesis had significantly higher 
detection rate (57%) at location 2 than those (7%) at loca¬ 


tion 10 (P = 0.0207) and Miniopterus fuliginosus showed 
significantly higher detection rate (50%) at location 10 than 
those in other locations (P = 0.0005), including 20% at 
location 1 and 0% at both location 5 and 7. 

Factors contributing for CoV detection were analysed and 
shown in Table 2. There was no sampling bias on genders 
because approximately half of the captured bats were females 
(49%) but the detection rate of CoV in females was 29%, sig¬ 
nificantly higher than that (17%) in males (P = 0.0348). 
Exceptionally, male Scotophilus kuhlii had a higher detection 
rate (50%) than females did (39%). By examining the nip¬ 
ples of female bats, 20 of 71 female bats evaluated have been 
lactating and four of them tested positive (20%) for CoV 
RNA. According to the descent of testicles into scrotum, 37 
of 67 male bats evaluated have been mating and only six of 
them were positive for CoV RNA (16%). 
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Table 2. Possible influence factors for coronavirus detection in Taiwan 


Factors 

Category 

% (CoV 

positive/total No.) 

P value 

Gender 

Female 

29% (35/121) 

0.0348 


Male 

17% (22/127) 


Age 

<1 year old 3 

8% (2/25) 

0.0778 


Adult 

25% (55/223) 


Lactation (Females) 

Lactating 

0% (0/6) 

0.1417 


Weaning 

29% (4/14) 



Non-lactating 

39% (20/51) 


Mating (Males) 

Mating 

19% (5/27) 

0.9137 


Prepare to mate 

10% (1/10) 



Non-mating 

20% (6/30) 



a Age of <1 year old was determined by epiphyseal-diaphyseal closure 
in the metacarpal-phalangeal joints of bats. 


Phylogenetic analysis of bat CoV sequences 

A 375-bp RdRp gene fragment from 24 bat CoV samples was 
sequenced and assigned GenBank accession number from 
KT381902 to KT381925 for bat CoV sample SI, S3, S5, S6, 
S8, S9, Sll, S13, S34, S37, S38, R6, R9, Rll, R14, M4, M8, 
Ml6, M22, M23, M24, M26, K9 and ML1 by the National 
Center for Biotechnology Information. The phylogenetic tree 
based on nucleotide sequences was constructed. Three dis¬ 
tinct lineages were defined and noted with the identities of 
nucleotide sequences to the closest CoV strains in Fig. 2. Bat 
CoV sample CYCU-R9/TW/2013 from Rhinolophus mono- 
ceros of Dongshan (location 10) belonged to SARSr CoV- 
related lineage in the genus Betacoronavirus because of the 
88% identity to those of human SARS Co Tor2 strain, bat 
CoV HKU3, and SARSr CoV Rp3 found in Rhinolophus bat 
of China. Scotophilus bat CoV 512-related lineage contained 
20 samples from four bat species and Miniopterus bat CoV 
lA-related lineage included three samples from two bat spe¬ 
cies. Both lineages belonged to the genus Alphacoronavirus. 
Bat CoV RNA isolated in all five Miniopterus bat samples 
from location 10 in southern Taiwan belonged to Scotophilus 
bat CoV 512-related lineage and those in two Miniopterus 
bat samples from location 2 in northern Taiwan belonged to 
Miniopterus bat CoV lA-related lineage. The same colony of 
Rhinolophus monoceros in the same irrigation culvert at loca¬ 
tion 10 was found bat CoV RNA belonging to SARSr CoV 
and Scotophilus bat CoV 512-related lineages. 

Discussion 

The current study has revealed a diverse bat CoV popula¬ 
tion circulating in Taiwan, following our previous report 
on Scotophilus bat CoV in central Taiwan (Su et al., 2016). 
Compared to the detection range of 5-73% for CoV RNA 
in bats worldwide, the CoV prevalence of 23% in Taiwan 
confirmed the widely distribution of CoV (Tang et al., 


2006; Tong et al., 2009; Rihtaric et al., 2010; Balboni 
et al., 2011; Shirato et al., 2012; Woo et al., 2012; Goes 
et al., 2013; Su et al., 2016). Female bats had significantly 
higher detection rate for bat CoV in this study, especially 
Rhinolophus monoceros (38% in females and 8% in males). 
In breeding season, female bats and their offspring tend to 
roost together as maternity colony and viruses are easier to 
circulate in large number (Gloza-Rausch et al., 2008). How¬ 
ever, male and female Scotophilus bats roost together so no 
significance was found in the detection rate of CoV RNA 
for different genders, reproductive and maturation status of 
Scotophilus bats in this study. The presence of bat CoV RNA 
in 17% of ground faecal samples could post a risk for poten¬ 
tial interspecies transmission through gene recombination, 
such as the emergence of human CoV (HCoV) NL63, 
HCoV HKU1, SARS CoV and MERS CoV (Bolles et al., 
2011), if the faeces containing bat CoV are ingested by ani¬ 
mals or humans. Further investigation on the survival of 
infectious bat CoV in the environment is needed. 

Both SARSr CoV and Scotophilus bat CoV 512 were 
found in the same colony of Rhinolophus monoceros. 
Although the bat detected positive for SARSr CoV 
appeared healthy, it had lower body weight of 4 g when the 
body weights of other collected Rhinolophus bats ranged 
from 4 to 6.5 g. Similar observation was noted in the Rhi¬ 
nolophus pearsoni carrying SARSr CoV Rp3 in Guangxi of 
China, which caused acute and self-limiting infection with 
viral clearance occurring between 2 weeks and 4 month 
(Lau et al., 2010). Recombination between SARSr CoV 
Rp3 from Guangxi and Rfl from Hubei may generate 
SARSr CoV SZ3 isolated from civets. A 4-year investigation 
on the migration pattern and genome epidemiology 
showed that Rhinolophus bats can migrate from 1.86 to 
17 km and carry bat CoV HKU2 or SARSr CoV Rp3 for 
further co-infection and frequent recombination (Lau 
et al., 2010). Scotophilus bat CoV 512 was only detected in 
one bat species, Scotophilus kuhlii , in Hainan Island of 
China and the Philippines previously (Tang et al., 2006; 
Watanabe et al., 2010). In this study, endemic Rhinolophus 
monoceros , endemic Kerivoula titania and migratory Min¬ 
iopterus fuliginosus were detected positive for bat CoV 
RNA belonging to Scotophilus bat CoV 512 related lineage 
beside Scotophilus kuhlii. Transmission of Scotophilus bat 
CoV 512 between migratory and endemic bat species was 
implied because some of them share the same roosting site. 

Bat CoVs sharing 96% of nucleotide identity to those of 
Miniopterus bat CoV 1A were not only detected in migra¬ 
tory Miniopterus fuliginosus but also endemic Submyotodon 
latirostris. It is also the first report that Miniopterus bat 
CoV 1 A can be detected in a non- Miniopterus bat species. 
These two bat species did not share roosting location 
although all bats tested positive came from northern 
Taiwan. Another unique finding is that two different bat 
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Fig. 2. Phylogenetic tree of three bat coronavirus lineages detected in Taiwan. The consensus 375 nucleotide sequences of partial RNA-dependent 
RNA polymerase gene from 24 bat CoV samples in Taiwan and selected sequences from Genbank were aligned using the Clustal W method and the 
tree was constructed using the neighbour-joining model with a 1000 bootstrap test and pairwise deletion option in the maximum likelihood method 
by mega6 program. Nucleotide identities were indicated next to the bracket. 


CoV lineages were detected in two different populations of 
Miniopterus fuliginosus whose living areas did not overlap. 
The distribution patterns of bat CoV in Miniopterus bats 
are very complicated. In Hong Kong, phylogenetically 


related Miniopterus bat CoV 1A and IB have distinct host 
specificities to Miniopterus magnate and Miniopterus 
pusillusy respectively, even though these two bat species co¬ 
inhabit the same cave. On the other hand, co-infection of 
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bat CoV IB and HKU8 was detected in Miniopterus pusil- 
lus. No evidence of recombination was detected between 
bat CoV 1A, IB, HKU8, 512, HKU2 and SARSr CoV iso¬ 
lated from Miniopterus bats (Chu et al., 2006). In Kenya, 
bat CoV 1A was isolated in Miniopterus africanus , Min¬ 
iopterus minor and Miniopterus natalensis , and both bat 
CoV 1A and HKU8 can be isolated in Miniopterus inflatus 
(Tong et al., 2009). In Bulgaria, bat CoV 1A, HKU-7, and 
HKU-8 were detected in Miniopterus schreibersii but not 
Rhinolophus and Hipposideros bats although all three bat 
species shared the same cave (Woo et al., 2012). Since Min¬ 
iopterus bats are migratory and the distance of migration is 
reported to be several hundreds of kilometers (Ramos Par- 
eira et al., 2009), it is implied that Miniopterus bat CoV 
might be transmitted between Hong Kong, Japan, even 
Europe as a result of bat migration. Further investigations 
on the distribution of bat CoV and interactions between 
migratory and endemic bat species are required for the 
close monitoring on SARSr-CoV and other bat CoV in Tai¬ 
wan. 

This study partially filled a knowledge gap regarding the 
prevalence and distribution of CoV in different bat species 
of Taiwan. Three bat CoV lineages identified as SARSr 
CoV, Scotophilus bat CoV 512 and Miniopterus bat CoV 1A, 
are very common in countries around Taiwan. Host switch 
may play a very important role in the distribution of differ¬ 
ent bat CoV from migratory bat species to Taiwan endemic 
bat species. Ongoing surveillance on bat CoV in different 
bat populations is necessary to understand the zoonotic 
potential of bat CoV. 

Acknowledgements 

This work was funded by the Ministry of Science and Tech¬ 
nology of Taiwan (NSC 102-2313-B-033-001). We thank 
Chia-Hong Chen at Shei-Pa National Park Headquarters at 
Miaoli for the foraging bat samples and Heng-Chia Chang 
for the ground faeces from Shuilin. 

References 

Agnihothram, S., R. Gopal, B. L. Yount, E. F. Donaldson, V. D. 
Menachery, R. L. Graham, T. D. Scobey, L. E. Gralinski, M. R. 
Denison, M. Zambon, and R. S. Baric, 2014: Evaluation of 
serologic and antigenic relationships between middle eastern 
respiratory syndrome coronavirus and other coronaviruses to 
develop vaccine platforms for the rapid response to emerging 
coronaviruses. /. Infect. Dis. 209, 995-1006. 

Annan, A., H. J. Baldwin, V. M. Corman, S. M. Klose, M. 

Owusu, E. E. Nkrumah, E. K. Badu, P. Anti, O. Agbenyega, B. 
Meyer, S. Oppong, Y. A. Sarkodie, E. K. Ealko, P. H. Lina, E. 
V. Godlevska, C. Reusken, A. Seebens, F. Gloza-Rausch, P. 
Vallo, M. Tschapka, C. Drosten, and J. F. Drexler, 2013: 


Human Betacoronavirus 2c EMC/2012-related viruses in bats, 
Ghana and Europe. Emerg. Infect. Dis. 19, 456-459. 

Balboni, A., A. Palladini, G. Bogliani, and M. Battilani, 2011: 
Detection of a virus related to Betacoronaviruses in Italian 
greater horseshoe bats. Epidemiol. Infect. 139, 216-219. 

Bennett, M., 2006: Bats and human emerging diseases. Epi¬ 
demiol. Infect. 134, 905-907. 

Bolles, M., E. Donaldson, and R. Baric, 2011: SARS-CoV and 
emergent coronaviruses: viral determinants of interspecies 
transmission. Curr. Opin. Virol. 1, 624-634. 

Casanova, L. M., S. Jeon, W. A. Rutala, D. J. Weber, and M. D. 
Sobsey, 2010: Effects of air temperature and relative humidity 
on coronavirus survival on surfaces. Appl. Environ. Microbiol. 
76, 2712-2717. 

Chu, D. K., L. L. Poon, K. H. Chan, H. Chen, Y. Guan, K. Y. 
Yuen, and J. S. Peiris, 2006: Coronaviruses in bent-winged 
bats ( Miniopterus spp.). /. Gen. Virol. 87, 2461-2466. 

Csorba, G., C. H. Chou, M. Ruedi, T. Gorfol, M. Motokawa, S. 
Wiantoro, V. D. Thong, N. T. Son, L. K. Lin, and N. Furey, 
2014: The reds and the yellows: a review of Asian Chrysop- 
teron Jentink, 1910 (Chiroptera: Vespertilionidae: Myotis). 

/. Mammal. 95, 663-678. 

Dobson, A. P., 2005: Virology. What links bats to emerging 
infectious diseases? Science 310, 628-629. 

Drexler, J. F., F. Gloza-Rausch, J. Glende, V. M. Corman, D. 
Muth, M. Goettsche, A. Seebens, M. Niedrig, S. Pfefferle, S. 
Yordanov, L. Zhelyazkov, U. Hermanns, P. Vallo, A. Luka- 
shev, M. A. Muller, H. Deng, G. Herrler, and C. Drosten, 

2010: Genomic characterization of severe acute respiratory 
syndrome-related coronavirus in European bats and classifica¬ 
tion of coronaviruses based on partial RNA-dependent RNA 
polymerase gene sequences. /. Virol. 84, 11336-11349. 

Ge, X. Y., J. L. Li, X. L. Yang, A. A. Chmura, G. Zhu, J. H. 
Epstein, J. K. Mazet, B. Hu, W. Zhang, C. Peng, Y. J. Zhang, 

C. M. Luo, B. Tan, N. Wang, Y. Zhu, G. Crameri, S. Y. Zhang, 
L. F. Wang, P. Daszak, and Z. L. Shi, 2013: Isolation and char¬ 
acterization of a bat SARS-like coronavirus that uses the 
ACE2 receptor. Nature 503, 535-538. 

Geller, C., M. Varbanov, and R. E. Duval, 2012: Human coron¬ 
aviruses: insights into environmental resistance and its influ¬ 
ence on the development of new antiseptic strategies. Viruses 
4, 3044-3068. 

Gloza-Rausch, F. A., A. Ipsen, M. Seebens, M. Gottsche, J. F. 
Panning, N. Drexler, A. Petersen, K. Annan, M. Grywna, S. 
Pfefferle. Muller, and C. Drosten, 2008: Detection and preva¬ 
lence patterns of group I coronaviruses in bats, northern Ger¬ 
many. Emerg. Infect. Dis. 14, 626-631. 

Goes, L. G., S. G. Ruvalcaba, A. A. Campos, L. H. Queiroz, C. de 
Carvalho, J. A. Jerez, E. L. Durigon, L. I. Davalos, and S. R. 
Dominguez, 2013: Novel bat coronaviruses, Brazil and Mex¬ 
ico. Emerg. Infect. Dis. 19, 1711-1713. 

Gouilh, M. A., S. J. Puechmaille, J. P. Gonzalez, E. Teeling, P. 
Kittayapong, and J. C. Manuguerra, 2011: SARS-Coronavirus 
ancestor’s foot-prints in South-East Asian bat colonies and 
the refuge theory. Infect. Genet. Evol. 11, 1690-1702. 


© 2016 Blackwell Verlag GmbH 


7 


Taiwan Bat Coronavirus 


Y.-N. Chen et al. 


Graham, R. L., and R. S. Baric, 2010: Recombination, reservoirs, 
and the modular spike: mechanisms of coronavirus cross¬ 
species transmission. /. Virol 84, 3134-3146. 

Ithete, N. L., S. Stoffberg, V. M. Corman, V. M. Cottontail, L. R. 
Richards, M. C. Schoeman, C. Drosten, J. F. Drexler, and W. 
Preiser, 2013: Close relative of human Middle East respiratory 
syndrome coronavirus in bat, South Africa. Emerg. Infect. Dis. 
19, 1697-1699. 

Ksiazek, T. G., D. Erdman, C. S. Goldsmith, S. R. Zaki, T. Peret, 
and S. Emery, 2003: A novel coronavirus associated with severe 
acute respiratory syndrome. N. Engl. J. Med. 348, 1953-1966. 

Lau, S. K., K. S. Li, Y. Huang, C. T. Shek, H. Tse, M. Wang, G. 

K. Choi, H. Xu, C. S. Lam, R. Guo, K. H. Chan, B. J. Zheng, 

P. C. Woo, and K. Y. Yuen, 2010: Ecoepidemiology and com¬ 
plete genome comparison of different strains of severe acute 
respiratory syndrome-related Rhinolophus bat coronavirus in 
China reveal bats as a reservoir for acute, self-limiting infec¬ 
tion that allows recombination events. /. Virol. 84, 2808-2819. 

Lau, S. K., K. S. Li, A. K. Tsang, C. S. Lam, S. Ahmed, H. Chen, 
K. H. Chan, P. C. Woo, and K. Y. Yuen, 2013: Genetic charac¬ 
terization of Betacoronavirus lineage c viruses in bats reveals 
marked sequence divergence in the spike protein of Pipistrel- 
lus bat coronavirus HKU5 in Japanese Pipistrelle: implications 
for the origin of the novel Middle East respiratory syndrome 
coronavirus. /. Virol. 87, 8638-8650. 

Lee, S. S., and N. S. Wong, 2015: Probable transmission chains 
of Middle East respiratory syndrome coronavirus and the 
multiple generations of secondary infection in South Korea. 
Int. J. Infect. Dis. 38, 65-67. 

Osborne, C., P. M. Cryan, T. J. O’Shea, L. M. Oko, C. Ndaluka, 
C. H. Calisher, A. D. Berglund, M. L. Klavetter, R. A. Bowen, 
K. V. Holmes, and S. R. Dominguez, 2011: Alphacoronaviruses 
in new world bats: prevalence, persistence, phylogeny, and 
potential for interaction with humans. PLoS One 6, el9156. 

Poon, L. L., D. K. Chu, K. H. Chan, O. K. Wong, T. M. Ellis, Y. 
H. Leung, S. K. Lau, P. C. Woo, K. Y. Suen, K. Y. Yuen, Y. 
Guan, and J. S. Peiris, 2005: Identification of a novel coron¬ 
avirus in bats. /. Virol. 79, 2001-2009. 

R Core Team (2013). R: A language and environment for statis¬ 
tical computing. R Foundation for Statistical Computing, 
Vienna, Austria. URL http://www.R_project.org/. 

Ramos Pareira, M. J., P. Salgueiro, L. Rodrigues, M. M. Coelho, 
and J. M. Palmeirim, 2009: Population structure of a cave¬ 
dwelling bat, Miniopterus schreibersii: does it reflect history 
and social organization? /. Hered. 100, 533-544. 

Rihtaric, D., P. Hostnik, A. Steyer, J. Grom, and I. Toplak, 2010: 
Identification of SARS-like coronaviruses in horseshoe bats ( Rhi¬ 
nolophus hipposideros) in Slovenia. Arch. Virol. 155, 507-514. 

Ruedi, M., G. Csorba, L. K. Lin, and C. H. Chou, 2015: Molecu¬ 
lar phylogeny and morphological revision of Myotis bats (Chi- 
roptera: Vespertilionidae) from Taiwan and adjacent China. 
Zootaxa 3920, 301-342. 

Shirato, K., K. Maeda, S. Tsuda, K. Suzuki, S. Watanabe, H. 
Shimoda, N. Ueda, K. Iha, S. Taniguchi, S. Kyuwa, D. Endoh, 
S. Matsuyama, I. Kurane, M. Saijo, S. Morikawa, Y. Yoshi- 


kawa, H. Akashi, and T. Mizutani, 2012: Detection of bat 
coronaviruses from Miniopterus fuliginosus in Japan. Virus 
Genes 44, 40-44. 

Su, B. G., H. C. Chen, H. C. Cheng, and Y. N. Chen, 2016: 
Detection of bat coronavirus and specific antibodies in 
chestnut bat ( Scotophilus kuhlii) population in central Taiwan. 
Taiwan Vet. J. 42, 1-8. 

Tang, X. C., J. X. Zhang, S. Y. Zhang, P. Wang, X. H. Fan, L. F. 
Li, G. Li, B. Q. Dong, W. Liu, C. L. Cheung, K. M. Xu, W. J. 
Song, D. Vijaykrishna, L. L. Poon, J. S. Peiris, G. J. Smith, H. 
Chen, and Y. Guan, 2006: Prevalence and genetic diversity of 
coronaviruses in bats from China. /. Virol. 80, 7481-7490. 
Tong, S., C. Conrardy, S. Ruone, I. V. Kuzmin, X. Guo, Y. Tao, 
M. Niezgoda, L. Haynes, B. Agwanda, R. F. Breiman, L. J. 
Anderson, and C. E. Rupprecht, 2009: Detection of novel 
SARS-like and other coronaviruses in bats from Kenya. Emerg. 
Infect. Dis. 15, 482-485. 

Watanabe, S., J. S. Masangkay, N. Nagata, S. Morikawa, T. 
Mizutani, S. Fukushi, P. Alviola, T. Omatsu, N. Ueda, K. Iha, 
S. Taniguchi, H. Fujii, S. Tsuda, M. Endoh, K. Kato, Y. Tohya, 
S. Kyuwa, Y. Yoshikawa, and H. Akashi, 2010: Bat coron¬ 
aviruses and experimental infection of bats, the Philippines. 
Emerg. Infect. Dis. 16, 1217-1223. 

Woo, P. C., S. K. Lau, C. S. Lam, C. C. Lau, A. K. Tsang, J. H. 
Lau, R. Bai, J. L. Teng, C. C. Tsang, M. Wang, B. J. Zheng, K. 
H. Chan, and K. Y. Yuen, 2012: Discovery of seven novel 
Mammalian and avian coronaviruses in the genus deltacoron- 
avirus supports bat coronaviruses as the gene source of alpha- 
coronavirus and betacoronavirus and avian coronaviruses as 
the gene source of gammacoronavirus and deltacoronavirus. 

J. Virol. 86, 3995-4008. 

Wu, Y., Y. Li, L. K. Lin, M. Harada, Z. Chen, and M. Motokawa, 
2012: New records of Kerivoula titania (Chiroptera: Vespertilion¬ 
idae) from Hainan Island and Taiwan. Mammal Study 37, 69-72. 
Yang, L., Z. Wu, X. Ren, F. Yang, J. Zhang, G. He, J. Dong, L. 
Sun, Y. Zhu, S. Zhang, and Q. Jin, 2014a: MERS-related Beta¬ 
coronavirus in Vespertilio superans bats, China. Emerg. Infect. 
Dis. 20, 1260-1262. 

Yang, Y., L. Du, C. Liu, L. Wang, C. Ma, J. Tang, R. S. Baric, S. 
Jiang, and F. Li, 2014b: Receptor usage and cell entry of bat 
coronavirus HKU4 provide insight into bat-to-human trans¬ 
mission of MERS coronavirus. Proc. Natl Acad. Sci. USA 111, 
12516-12521. 

Zaki, A. M., S. van Boheemen, T. M. Bestebroer, A. D. Oster- 
haus, and R. A. Fouchier, 2012: Isolation of a novel coron¬ 
avirus from a man with pneumonia in Saudi Arabia. N. Engl. 

J. Med. 367, 1814-1820. 

Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Table SI. The 28 coronavirus sequences from NCBI for 
phylogenetic analysis. 


8 


©2016 Blackwell Verlag GmbH 


